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We study numerically the non-equilibrium critical properties of the Ising model defined on direct 
products of graphs, obtained from factor graphs without phase transition (Tc = 0). On this class of 
product graphs, the Ising model features a finite temperature phase transition, and we find a pattern 
of scaling behaviors analogous to the one known on regular lattices: Observables take a scaling form 
in terms of a function L(t) of time, with the meaning of a growing length inside which a coherent 
fractal structure, the critical state, is progressively formed. Computing universal quantities, such as 
the critical exponents and the limiting fluctuation-dissipation ratio Xao, allows us to comment on 
the possibility to extend universality concepts to the critical behavior on inhomogeneous substrates. 
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I. INTRODUCTION 

The equilibrium physics of second order phase transitions is quite well understood nowadays, due primarily to the 
development of scaling theories and of the renormalization group. Systems above the lower critical dimension build 
up an increasing coherence length ^ as a finite temperature is approached. In the neighborhood of Tc, ^ has 
grown much larger than any other characteristic length, and physical quantities can be expressed in scaling forms 
in terms of ^ only. Universal quantities, such as the exponents entering those relations, are known to depend only 
on a small set of parameters as, for instance, the space dimensionality d and the symmetry of the order parameter. 
Other features, among which the local geometry of the underlying lattice (e.g. triangular, cubic etc ...), are known to 
affect only non-universal quantities, e.g. the value of Tc. Scaling and universality concepts can be extended also to 
far-from-cquilibrium systems approaching the critical state kinetically, as in the prototypical case of samples quenched 
from high temperature to the critical temperature T^. In this case an infinite coherence length is built out of a length 
L{t) which grows indefinitely in time, and scaling, with respect to L{t) in this case, is again observed. 

This general picture of critical phenomena is well established when the underlying lattice is homogeneous, and 
interesting results have been obtained in the case of fractal structures, exploiting their scale invariance and their 
properties under exact decimation procedures However a general discrete structures, i.e. a graph, which can feature 
strong inhomogeneities and have no a priori symmetry, is defined by a purely topological information encoded in the 
nearest neighborhood relations between its sites, and in this case the situation is still unclear. This happens despite 
examples of inhomogeneous graphs, ranging from disordered materials, to percolation clusters, glasses, polymers, and 
bio-molecules, may be abundantly found in physics, economics, chemistry and biology 0. 

A first question concerns the topological properties playing the role of the Euclidean dimension d in determining 
the existence of a (finite temperature) phase transition. In the case of continuous symmetry models [1, Ql (e.g. the 
Heisenberg model) this role has been shown to be played by the spectral dimension dg of the graph, a quantity related 
to the low eigenvalues behavior of the density of states in the Laplacian operator Q . The spectral dimension of a 
graph alone determines the existence of a finite temperature phase transition, and ds = d^ = 2 i& the lower critical 
dimension above which a graph can sustain an ordered phase at finite temperature T < Tc. Therefore one expects 
a pattern of scaling behaviors similar to the one observed on homogeneous lattices and a natural extension of the 
concept of universality, with ds playing a role analogous to the Euclidean dimension d. However, the relation between 
the spectral dimension and the critical exponents in the case of continuous symmetry models on a general graph is 
still an open problem, with a few rigorous results 0, Q . 

For discrete symmetry models, even such a simple topological generalization of the Euclidean dimension does not 
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exists, and an indicator playing the same role as dg in this case is not known, although partial results have been 
obtained in this direction ■ Discrete symmetry models have been shown to feature a phase transition when dg > 2 
[3] but a necessary and sufficient condition for its existence has not been yet demonstrated, and also the effect of 
topology on critical exponents is still a completely open problem. In particular, the knowledge of 0?^ is not sufficient 
to determine the universality class of a discrete symmetry model. This poses a number of important questions yet 
to be answered about the critical properties of discrete symmetry models on general graphs. Lacking an equivalent 
of ds, the problem of distinguishing graphs where a phase-transition occurs at finite Tc from those where Tc = 0, is 
opened. 

Let us recall that on usual lattices the lower critical dimension of the model is = 1, whereas a two-dimensional 
lattice with Tc > can be obtained as the direct product (as defined in Sec. |IT]) of two one-dimensional systems. 
Extending this observation to the realm of inhomogeneous structure and general graphs, one is led to consider direct 
products of graphs with Tc = 0. 

The direct product of graphs represents a practical receipe to build a class of inhomogeneous structures with Tc > 
but no a priori symmetry, allowing the study of their critical properties as a function of topology alone. Following 
this reasoning, in this paper we study numerically the non equilibrium critical properties of the Ising model defined 
on direct products graphs. We determine the numerical value of Tc, which is finite as predicted by the analytic results 
Q and we show that the whole non-equilibrium scaling behavior is analogous to that found on usual lattices with 
T > Tc- In particular, observables as the two-site/two-time correlation function Cij{t, s) can be expressed in terms of 
a growing length L(t) inside which a coherent fractal structure, the critical state, is progressively built. 

The behavior of the Ising model on fractal structures with Tc = 0, considered previously in [8|, is very similar to 
that observed on the usual 1 — d lattice. Not only on these fractal structures the model remains disordered at any 
finite temperature but, in addition, at least a couple of universal quantities (the response function exponent a of Eq. 
, and the limiting fluctuation-dissipation ratio Xoo of Eq. ([T7)) ') take always the same value of the 1 — d lattice 
This features might be interpreted as an indication of a sort of universality, although in a vague and weak sense, 
since other exponents are different. This makes the direct products of these Id-like graphs interesting also because, 
pushing the above presumptive universality arguments even further, and in analogy to what happens on homogeneous 
lattices, one might wonder if the 2d-like graphs obtained in this way show any universal behavior, at least in the weak 
sense discussed above. However, the numerical results of this paper indicate that all the scaling exponents, including 
a, and also X^c, differ among the possible product graphs and with respect to the values taken in the homogeneous 
case with d — 2. This proves that some other topological differences are relevant (in a renormalization group sense), 
making the concept of universality, if any, yet obscure. Understanding the nature of these relevant features, and, 
possibly, the way to associate them some topological indexes, analogous to d and dg, remains an open problem. 

The paper is organized as follows: In Sec. |lT] we recall some basic notions about graphs. Sec. IIIII contains an 
overview of the scaling behavior on homogeneous structures, and extensions to the realm of the inhomogeneous ones. 
In Sec. IIVI we study numerically the critical properties of the Ising model on some product graphs, determining the 
critical temperature and the scaling properties of correlation and response functions. Sec. |V] contains a brief analyses 
of the kinetics of the model after a quench below Tc, in order to comment on the issue of universality also in this case. 
In Sec. IVII we draw our conclusions and discuss some open problems. 




II. GENERALITIES ON GRAPHS 



A graph (network) Q is defined by a countable set of sites i connected pairwise by unoriented links The 
chemical distance r^j Q, i.e. the number of links of the shortest path connecting sites i and j, defines a natural 
metric on Q. The van Hove sphere So,r of radius r and center o in this metric is the sub-graph of Q composed by the 
sites whose distance from o is smaller than r + 1. We call No^r the number of sites in So.r- On infinite graphs the 
asymptotic behavior for large r of Ng^r defines the fractal dimension: 

No^r r'^f (1) 

where ^ denotes the behavior for large r. In the following we will consider only connected graphs embeddable in a 
finite dimensional space, with df well defined and finite. We also require that the degree Zi (number of neighbors of 
the site i) is bounded. On fractal graphs one can cquivalently explore the infinite structure using, in a very natural 
way, finite generation (sub)fractals instead of the Van Hove spheres (see below for a definition of generation). 

The graph can be described using a set of characteristic matrices. The adjacency matrix Aij of a graph has entries 
equal to 1 if i and j are neighboring sites ({«,.?} is a link) and Ai_j = otherwise. The Laplacian matrix A^j- is 
defined as 

^i.j — ^i.j^i -^i-j (2) 
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where Zi = Ai_j is the degree of i. Interestingly, Aij is the generahzation to graphs of the usual Laplacian 
operator of Euclidean structures [T]| . In particular its spectrum is positive and for connected structures the constant 
vector is the only eigenvector of eigenvalue zero. We define the spectral density of the subgraph So,r as Pr{l) = 

N~l ^1,^=1 ~ Lk,r) where L^^r are the eigenvalues of the Laplacian matrix of So.r- Let p{l) be the limit of Pr{l) 
for r — )■ oo . For positive functions, if the thermodynamic limit r oo exists, then it can be shown that the result is 
independent of the center of the sphere o If p{l) behaves for small / as 

p{l) ^ (3) 

then ds is defined as the spectral dimension of the graph Q . 

Graphs of large dimension can be built by introducing the direct product. Given two graphs Q and T-L, the direct 
product X H is a graph whose sites are labeled by a pair with i and j belonging to Q and % respectively, 

(i, j) and are neighbor sites in Q x H if i ~ i' and {j, j'} is a link of H, or if j = j' and {i,i'} is a link of 

Q. Interestingly, a basic properties oi Q x H is that, calling and d^ the dimensions of the graph Q, the following 
relations hold [ill 

rff^^-df+dl' df^^d^^+dj (4) 

i.e. the spectral and fractal dimensions of the product graph are the sum of the dimensions of the original graphs. 

We will consider models defined on graphs with known d^ and d/ to verify the relevance of these dimensions on 
universal properties and we will use direct products to obtain graphs with a finite critical temperature. In particular, 
we use as base graphs two finitely ramified fractals [Hj with ds < 2, the T- fractal (TF) and the Sierpinski gasket 
(SG). These fractals can be built recursively as shown in Fig. [TJ For the SG one starts with three spins located on 
the vertexes of a triangle. This is the first generation Gi of the structure. The second generation G2 is obtained by 
considering three adjacent Gi structures, building a larger triangular object. The process is then iterated at will. A 
similar construction is used for the TF, as shown in Fig. [l] the first generation is Gi and the generation Gn can be 
built attaching 3 G„_i structures in a single site. We will denote with i?„ the size of the finite generation G„ i.e. 
the maximum distance between two sites of G„. The fractal and spectral dimensions of the TF and of the SG can 



be analytically evaluated by means of exact renormalizations [l3[, yielding df = log(3)/ log(2), dg ~ log(9)/ log(6) 
and df ~ log(3)/ log(2), dg = log(9)/ log(5) respectively. We consider the graphs obtained from the product of two 
SG's (SGxSG) and two TF's (TFxTF), their fractal and spectral dimensions can be calculated from Equations @ 
obtaining df = 3.17.., 4 = 2.45... for the TFxTF and df = 3.17.., 4 = 2.73... for the SGxSG. Since ds > 2 the Ising 
model present a phase transition at finite temperature Tc, according to the generalized Froelich-Simon-Spencer bound 
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III. MODELS AND SCALING RELATION 



In this section we introduce the model, the dynamical quantities of interest and overview their known scaling 
behavior. We consider, on a given graph Q, the Ising model, defined by the Hamiltonian 

H[a] = - J o-i<^j = - J ^^Ajcr^aj (5) 

<ij> i,j 

where ct,; = ±1 denotes the spin of site i, < ij > are nearest neighbors on the graph and Aij is the adjacency matrix. 
The dynamics is introduced by randomly choosing a single spin and updating it with Metropolis transition rates: 

w{[c7] ^ [a']) = min [1 , cxpi- AE / {K bT))] . (6) 

Here [a] and [a'] are the spin configurations before and after the move, Kb is the Boltzmann constant, and 

AE = H[a'] - H[c7]. (7) 

A Montccarlo time step is a sequence of random moves where N is the number of sites in the finite realization 
of the graph Q. We will focus on the non-equilibrium situation where a completely random initial configuration, 
corresponding to an equilibrium state at T = oo, is instantaneously quenched to a temperature equal or smaller than 



A. Quenches to Tc 



After quenching our system to the critical temperature, the two-time correlation function in the case of Ising 
variables is defined as 



Cij{t,s) = {ai{t)aj{s)) 



(8) 



((. . .) being an average over initial conditions and thermal histories). On a homogeneous structure (i.e. a lattice) for 
sufficiently large times ^ takes the scaling form [13, [HI 



Cy(t,s) = L(s)-''^5 



HsYLis) 



(9) 



where L{t) oc t^^^ is the typical length associated to the growth of the critical phase, z is the usual equilibrium 
dynamic exponent relating the relaxation time to the coherence length r oc r is the Euclidean distance between 
sites i and j, g{x,y) is a scaling function, and 6z = d — 2 -I- is an exponent related to the usual equilibrium critical 
one rj. The scaling form (j9|) is determined by the growth of correlated regions of size L{t) with a fractal dimension 
Df = d — zb/2 [l^. Notice that when d — >■ d^, Tc ^ and 6 — > 0. Then Df — > d and these regions grow compact. 
For instance, for the Ising model we are considering here one has 77 = 1 at d = = 1, leading to 6 = 0. This is a 
manifestation of the fact that the quench at is not a critical quench, but rather resembles a quench in the ordered 
phase, as it will be discussed in Sec. IIIIBI Let us remark that the analysis of the out of equilibrium process provides 
the equilibrium static and dynamic exponents, since they enter the form ([9]). 

In the case of inhomogeneous structures, one might still expect some sort of scaling to hold. However, in that case 
a unique definition of distance r is not available. Hence, in what follows, it will be useful to introduce the space 
integrated correlation function 



(10) 



i.'jeG 



where G is a certain subset of sites, which has the advantage of a straightforward generalization to fractal structures. 
Notice that we have restricted the definition to the equal time correlation s ~ t, since only this function will be 
considered in Sec. IIVI In the case of homogeneous lattices, considering a box G of size R, using Eq. ((9)) and 
L{t) cx t^^^, one finds 



Fait) = t 



-b+d/z 



f 



R 
-jTTz 



(11) 
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On an inhomogcneous structure such as the TFxTF or the SGxSG we will consider the quantity 

FG,{t)= J2 '^^M^ (12) 

where sites are summed over the internal sites of the fc-th generation Gk , for which we expect the scaling 



(13) 



On fractals the scaling hypothesis is hence naturally defined by Eq. (fT3|) where the scaling function / depends on 
Rk/t^^^ only. Notice that the use of Eq. ([T3|) avoids the problem of the definition of a distance, since the length i?, 



fe 



is naturally associated to the fc-th generation. 

In Sec. pvp we will study in detail the so called autocorrelation function, obtained by letting i = j in Eq. 
We do this not only because it is very often considered in aging systems, but also because, being an on-site quantity, 
it circumvents the definition of distances. In homogeneous systems this quantity does not depend on position and, 
denoting it as C(t, s), from Eq. ([9]) and L{t) oc one has the scaling behavior 



C{t,s)^s~''h[-\, (14) 



with h{t/s) ~ g[0, L{t) / L{s)]. In equilibrium conditions, C{t,s) is associated to the autoresponse function R{t,s) = 
5{(Ji) / 5hi{s)\f^.^Q, describing the effect of an impulsive perturbing magnetic field hi{s) switched on at time s on site i, 
by the fluctuation-dissipation theorem TR{t, s) = —dC{t, s)/ds. When the system is out of equilibrium after a quench 
to Tc, the fluctuation-dissipation theorem no longer holds in general. However, in the short time difference regime, 
t/s ~ 1, due to local equilibrium, observables such as C and R behave as in equilibrium [TtI . I18| . The constraint 
imposed by the fluctuation-dissipation theorem, then, determines the value a = 6 of the exponent entering the scaling 
form [ll [ll [13 d of i? 

R{t,s)^s-^''+^^h(^^y (15) 
With the scalings (|14I15|) and a = b, the fiuctuation-dissipation ratio 

^^(^)-Tm (16) 



is a function of t/s. Its limiting value 



dC(t 
ds 



Xoo = lim lim X (-] (17) 



is of a particular interest, since it was shown [IJ, 1T5| to be an universal quantity. For the Ising model in d = 2 one has 
~ 0.33 p^ . On inhomogcneous structures, since the autocorrelation function and the autoresponse may be site 
dependent, we define them as C{t, s) — (1/^) J2i Cii{t, s) and similarly for R{t, s), where N is the number of sites of 
the finite realization of the infinite graph Q. 



B. Quenches below Tc 



When the quench is performed below the nature of the process changes, because the target equilibrium state is 
ordered (magnetized) and degenerate. In this case domains of the two possible equilibrium phases grow, and their 
geometry is compact. The observable quantities introduced above split into two contributions [2Q] 

Cit,s)^Cst{t-s) + Cagit,s), (18) 

and similarly for the other ones. The first stationary term is a quasi-equilibrium contribution provided by the interior 
of the domains, which is in fact equilibrated in one of the two possible phases. What is left over, the non-equilibrium 
character due to presence of the interfaces, contributes with the second term, which takes the scaling form (jl4p with 
b = (due to the compact geometry of the domains and the above-mentioned relation between b and Df). Similarly, 
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for the aging part of the response function one has scahng as in Eq. ((T5|) . At variance with critical quenches, there is 
no known relation between the exponent a and 6, nor with any other equilibrium or non-equilibrium exponent. In two 
dimensions the value a — 1/4 has been conjectured, and numerical simulations tend to conform to this hypotheses 
PH (although the value a = 1/2 has been also reported [13 )■ It is also known that a > for d > di,, and that a ^ 
when d ^ cLl. Hence, the value of the exponents a and h, when d is lowered towards d^ approach the value a = = 
both along the path of critical quenches (at T = Tc) or along the route with T = (or any route with T < Tc). 
However, as explained in [l7| the quench ai d — dL-^T = {) cannot be regarded as critical, but rather as a quench in 
the ordered phase, as it is clear since the equilibrium state at T = is degenerate and magnetized. Finally, let us 
recall that a > and 6 = implies that X^o = through Eqs. (|16ll7p . 

C. What is known about a and Xoo on inhomogeneous structures 

We have already discussed the fact that on Euclidean geometries one always finds a > in the phase-ordering kinetics 
of systems above d^, while a = at d^ [l3, H^. Elaborating on this, in d, [gl] it was conjectured that a similar 
property holds also on inhomogeneous structures, and the exponent a was proposed to infer the presence/absence of 
a finite temperature transition on a graph. Indeed a = was found on structures such as the SG, the TF and others 
where a magnetized state cannot be sustained at any finite temperature, whereas a > was observed for instance 
on the toblerone lattice (the product graph between the SG and the line) where Tc > 0. Regarding the limiting 
fluctuation-dissipation ratio, the data reported in yield results very well compatible with the value Xoo = 1/2 
observed in the d = 1 homogeneous case [25| . on all the structures with a = considered. Apart from a and X^o, 
other exponents such as z were found to be temperature dependent. In conclusion, all the graphs with Tc = 
considered in share the same value for two universal quantities, a and Xoo, with the usual Id lattice. 

IV. NUMERICAL RESULTS: QUENCH AT n 

In the following we will study the critical behavior of the product graphs TFxTF and SGxSG. In the simulations 
we used finite fractals of generation 7 so that the total number of sites equals 4787344 for the TFxTF and 1199025 
for the SGxSG. We verified that the graph size is large enough to avoid finite size effects so that average quantities 
such as C{t, s) turn out to be independent of the size N. In the following we will set Kb — 1 and J = 1. We will 
consider first the scaling of the equal time correlation FGf^{t), which, as a byproduct, allows one to determine Tc, and 
then the properties of two-time quantities. 

A. Determination of Tc and scaling of Fq^. (t) 

In order to determine Tc we have followed two methods. The first makes use of the critical finite-size scaling of 
Fg^ (<), Eq. ([T^]). This form holds right at Tc, with 6 > 0. Instead, above Tc there is a departure from scaling at large 
times, since the system eventually equilibrates. Below Tc, a similar scaling form holds but only for the aging part of 
the correlation (recalling the discussion around Eq. ^TE\\ ). and with the value b = characteristic of phase-ordering. 
In conclusion one expects to observe deviations from the scaling (|13p with 6 > as one moves away from Tc- Hence, 
the method we use for determining Tc amounts to the computation of Fcf, (t) at several temperatures, trying to 
collapse the data according to Eq. by plotting ''■>'/ FQf^{t) against Rk/t^^^, using T, & > and z are fitting 
parameters. More precisely, using the least square method, we introduce a quantity S{T, z, b) with the meaning of a 
variance around the optimal value (see the Appendix for a precise definition). S{T, z, b) measures the quality of the 
collapse and then we determine Tc, z and b as the values that provide the minimum of S{T, z,b). Examples of data 
collapsed in this way arc shown in Fig. [2] for the TFxTF. In the three panels we show the behavior of the curves 
collapsed with the optimal choice of Tc, z, b (upper panel) or with other choices of T (larger and smaller than the 
optimal one, which turns out to be Tc — 3.21) for which S{T, z, b) is slightly larger than the minimum. It is clearly 
observed that the collapse is better in the first case. 

The next point is to give an estimate of the errors. It must be stressed here that statistical errors are quite small in 
our simulations, whereas the major source of inaccuracy is a systematic effect due to having a finite window of times t 
and of sizes Rk- Therefore, in order to estimate the errors on the fitting parameters Tc, z, b, we define a compatibility 
region, namely a region around the optimal value of the parameters where the value of S{T, z, b) is not larger than 
n times the minimum. Here n is the tolerance that we fix to n = 3 (this value was chosen arbitrarily, but we have 
checked that different choices, ranging from ri = 2 to n = 5, do not change the compatibility region significatively). 
An example of this procedure is depicted in Fig. [3] (left panel). Here, since we have three parameters, for the purpose 
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FIG. 2: (Color online.) Data collapse obtained by plotting t'' ''f^^Fcit) against Rk/t^''^ for the TFxTF at three temperatures, 
T — 3.21 ~ Tc (upper panel), T = 3.19 (lower left panel) and T = 3.23 (lower right panel). Different points are evaluated on 
TFxTF of generation k with k ranging from 3 to 7. The value of S(T, z, &) is 6 ■ 10"'' at T = 3.21, 3 ■ 10"^ at T = 3.19 and 
2 ■ 10"^ at T = 3.23 
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FIG. 3: (Color online.) In the left panel, continuous lines enclose the parameter compatibility regions (see text) for the TFxTF 
and the SGxSG. The right panel represents the data collapse for Fof, (t) on the SGXSG. Upper plot is obtained by using as 
scaling exponents the value of the center of the compatibility region, evidencing a very good data collapse. The lower panel 
refers to values of the exponents at the border of the compatibility region. Notice that already for n = 3 the scaling is sensitively 
worse, this evidence the reliability of the definition of the compatibility region. 



of visualization wc have fixed one of them, the temperature, to its optimal value and we have plotted the projection 
of the compatibility regions (for the TFxTF and SGxSG) on the z, b plane. The right panels show that the collapse of 
the curves for the SGxSG gets worst when one moves from the optimal value of the parameters, in the centre of the 
compatibility region, to the border. From the determination of the compatibility region we can eventually associate 
an error to the determination of T^, z, 6 as the distance between the optimal value and the border of the region. This 
whole procedure provides the following values = 3.21 ± 0.02, 6 = 0.4 ± 0.15, z = 2.6 ± 0.6 for the TFxTF, and 
Tc = 6.14 ± 0.02, 6 = 0.3 ± 0.15, z = 2.8 ± 0.6 for the SGxSG. The errors on z and b are quite large, therefore within 
this approach it is not clear if the critical behavior is characterized by the same exponents for our different structures. 
However Fig. [3] evidences that it is very unlikely that both exponents are equal because the compatibility regions do 
not intersect. 

Comparing the values of Tc for the TFxTF and the SGxSG with the values of the lattice in two, three and 
four dimensions (i.e. Tc ~ 2.27, Tc ~ 4.51 and Tc ~ 6.68) one observes that Tc tends to increase with the average 
coordination number of the graph. However, clearly other non-universal parameters play a relevant role in determining 
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FIG. 4: (Color online.) Left panel: Plot of C{t,s) against s for different fixed values of t/s for the TFxTF quenched to Tc- 
The continuous line is the power law with b = 0.33. Right panel: Data collapse obtained by plotting s^C{t, s) against t/s, 
for several values of s. 



the critical temperature, as can be noticed from the fact that the SGxSG and the 4-dimensional lattice have both 
Zi ~ % but a different value of T^- 

The second method for the determination of Tc makes use of the scaling properties of C{t,s), as proposed and 
discussed in [2^. With this technique we obtained results in agreement with the first method. Let us stress the 
advantage of these non-equilibrium methods, since there is no need to equilibrate the system which, due to the critical 
slowing down, is a very demanding numerical task. 



B. Scaling of two-time quantities 

Letting T = Tc, we have computed the two-time quantities C{t,s) and R{t,s). In view of Eq. (|14p. by plotting 
C{t, s) for fixed values of t/s against s, (left panel of Fig. |3]and Fig. [5]) we have determined the value h = 0.33±0.03 for 
the TFxTF and h = 0.39 ±0.03 for the SGxSG, in good agreement with the previous determination obtained through 
Fc{t). Notice that the two time approach provides a much more precise value of b: indeed a single scaling exponent 
have to be fitted while in the scaling approach to Fckii) both the exponents b and 1/z are evaluated numerically. 
Reconsidering the results of Sec. IIV Al we obtain also a better estimate of the z's, which turn out to be z = 2.9 ± 0.3 
and z = 2.4 ± 0.2 for the TFxTF and SGxSG respectively. Notice that these values (particularly b) are very different 
from those known for the homogeneous 2d lattice, namely b ~ 0.115 and z ~ 2.17. With a good confidence one may 
also conclude that they are different also for the two product graphs (although, in principle, the values b = 0.36, 
z = 2.2, obtained at the upper/lower limit of the error bars could be compatible with both the structures). The data 
collapse obtained by plotting s^C{t, s) against t/s can be checked in the right panel of Fig. 2] and [51 For large values 
of t/s the collapse is very good. Deviations are observed in the small t/s- region, but these appear to become less 
important as s increases. This is a quite common feature in phase-ordering, observed also on homogeneous lattices, 
and can be attributed to pre-asymptotic corrections [3 HH [13 ■ 

For the computation of the response function we have used the field- free algorithm introduced in [2^. We have 
checked that the scaling form ([15]) is obeyed with an exponent consistent with the expected behavior a = b. In Fig. 
[HI we plot the fluctuation-dissipation ratio against s/t, for different values of s. The good data collapse observed 
confirms that R scales as in Eq. (|15p . with a ~ b. For s/t :^ 1, in the quasi-equilibrium regime, there are no sensible 
deviations from the fluctuation-dissipation theorem, as discussed in Sec. IIII Al and X ~ 1. For s/t <C 1 the so called 
aging regime is accessed, ad X lowers. Extrapolating the intercept at s/t = we obtain X^o ~ 0.36 ± 0.03 for the 
TFxTF and X^o = 0.40 ±0.03 for the SGxSG. Notice that these values are compatible within error bars. However, by 
comparing these values with the one {Xoo — 0.33) measured in the usual square lattice, one sees that for the TFxTF 
they are only marginally compatible, while for the SGxSG they are totally incompatible. 

In conclusion, in the case of critical quenches of the Ising model on the product structures considered here one flnds 
a finite Tc and a pattern of behaviors analogous to that found on homogeneous lattices. The value of the universal 
(on usual lattices) quantities a and X^o, appear to be different between the TFxTF, the SGxSG and the 2d lattice. 
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iog,o(s) iog,(,(t/s) 

FIG. 5: (Color online.) Left panel: Plot of C{t, s) against s for different fixed values of t/s for the SGxSG quenched to Tc. The 
continuous line is the power law with b = 0.39. Right panel: Data collapse obtained by plotting s''C{t,s) against t/s, for 
several values of s. 
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FIG. 6: (Color onhne.) Plot of X{t,s) against s/t for different fixed values of s for the TFxTF (left) and SGXSG (right) 
quenched to Tc. 

V. NUMERICAL RESULTS: QUENCH BELOW Tc 

For completeness, for the TFxTF we have computed the value of the exponent a also in the case of a quench 
below Tc (XoD is trivially zero in this process, see discussion at the end of Sec. IIII B[) . It must be recalled that in 
the Toblerone lattice an exponent a was found, whose value was compatible with the value a = 1/4 found in the 
two-dimensional case. This could suggest that the conjecture discussed above, namely that products of Id-like graphs 
could form a sort of universality class sharing the exponent a, although not confirmed by our previous data in critical 
quenches, could hold at least restricting to subcritical quenches. In order to check this issue we have computed the 
integrated response function x(f, s) = J* dt'R{t, t') and, recalling the additive form discussed in Sec. IIIIBl we have 
isolated the aging term Xag by using a dynamics where flips of spin in the bulk arc forbidden, as discussed in [2]| . 
From Eq. ([TS]) one has 

Xaa{t,s)^s~''h{t/s), (19) 

where h{x) is another scaling function. The data of Fig. ^ show that a good scaling collapse is obtained with an 
exponent a ~ 0.13 ± .02 which is different from the value found in d = 2 (see inset of Fig. ([7]) ). The fact that a 
takes a comparable value in the = 2 homogeneous lattice and on the Toblerone lattice, therefore, seems not to be a 
general property of product graphs. 
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FIG. 7: (Color online.) s^xC^jS)! with a = 0.13, is plotted against t/s for different fixed values of s for the TFxTF quenched 
to T = 3.0 < Tc- In the inset we plot s^'^^xC^) against t/s showing that for a — 1/4 the scaling is not satisfied. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper we have studied the scaling properties of the Ising model quenched to or below the critical temperature 
on graphs obtained by making direct products of Id-like structures, such as the TFxTF and the SGxSG. The direct 
product is a convenient tool to build graphs with a topology sustaining a finite Tc, where the critical properties can be 
studied. This allows one to investigate the critical properties of such structures and to determine the critical exponents. 
Moreover, dynamical aspects can also be considered, by studying the evolution after temperature quenches. The aim 
of our analysis is to check the validity of dynamical scaling on inhomogeneous structures, and to study the behavior of 
the (limiting) fluctuation-dissipation ratios and of the scaling exponents. On regular lattices these are well understood 
and their universal properties are well known. On product graphs we found a non-equilibrium scaling behavior similar 
to that found on homogeneous lattices above d^, where time enters observable quantities through a single growing 
length representing the size of the critical correlations established. Regarding the quantities which on homogeneous 
lattices are universal, previous studies of phase-ordering on a certain class of graphs showed their dependence on 
several parameters, among which the temperature, at variance with the regularity observed on homogeneous lattices. 
A notable exception was represented by the exponent a and the limiting fluctuation-dissipation ratio X^o , which were 
found to take the same values a = and X^o = 1/2 on all the graphs with Tc — Q considered insofar. In this paper 
we have considered the possibility that such a regularity could be extended to the various direct products of these 
Tc = graphs than one can consider, by studying if critical exponents or Xoo take a unique value for all the graphs 
of this class. The results we found, however, are negative in this respect, either in critical quenches or in sub-critical 
quenches. Indeed we have shown that product graphs exhibit different exponents, including a, and, restricting to 
critical quenches, also a different value of X^o- This may indicate that a robust universality property as observed 
on usual lattices is lost in the realm of inhomogeneous structure, or that, if some analogue of it exists, the relevant 
parameters of the graph topology determining universality have not yet been identified. 



Acknowledgments 



F.Corberi acknowledges financial support from PRIN 2007 JHLPEZ [Statistical Physics of Strongly Correlated 
Systems in Equilibrium and Out of Equilibrium: Exact Results and Field Theory Methods). 



Appendix 



We define the quality of the collapse between the different curves as follows. First, for every time ti we fit the curve 
log(FG^ [ti)) as a function of the size Rk with a polynomial of the form F{ti, x) where x = \og{Rk/t- ). This is done 

since the function Fc^(t) is known only on a discrete set of sizes Rk- The curves F[x,ti) = \og{t\ '^^1^) \- F{ti,x), 
obtained at different times ti should collapse once plotted against x (the reason for taking the logarithms is that. 
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since one has power-law dependences, the logarithms allow one to better take into account deviations from perfect 
collapse on a wide range of sizes and times). Then we introduce a squared distance between the rescaled curves as 
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